Anomalous and dimensional scaling in anisotropic turbulence 
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We present a numerical study of anisotropic statistical fluctuations in homogeneous turbulent 
flows. We give an argument to predict the dimensional scaling exponents, Cd(p) — (p + i)/3, for 
the projections oi p-th order structure function in the j-th sector of the rotational group. We show 
that measured exponents are anomalous, showing a clear deviation from the dimensional prediction. 
Dimensional scaling is subleading and it is recovered only after a random reshuffling of all velocity 
phases, in the stationary ensemble. This supports the idea that anomalous scaling is the result of a 
genuine inertial evolution, independent of large-scale behavior. 



In recent years a huge amount of theoretical, numer- 
ical and experimental work has been done in order to 
study anisotropic turbulent fluctuations Typical 
questions go from the theoretical point of calculating and 
measuring anomalous scaling exponents in anisotropic 
sectors |4[-[7[|l]|, to the more applied problem of quan- 
tifying the rate of recovery of isotropy at scales small 
enough Another important issue is the univer- 

sality of anisotropic scaling exponents, i.e. whether they 
are an intrinsic characteristic of the Navier-Stokes non- 
linear evolution or they are fixed by a dimensional match- 
ing with the external anisotropic forcing. 

Important steps forward in the analysis of anisotropic 
fluctuations have recently been done in Kraichnan mod- 
els, i.e. passive scalars/ vectors advected by isotropic, 
Gaussian and white- in-time velocity fields |l3[ , with a 
large-scale anisotropic forcing ||l4-16|. In those models, 
anomalous scaling arises as the result of a non-trivial null- 
space structure for the advecting operator (zero modes). 
Also, correlation functions in different sectors of the rota- 
tional group show different scaling properties. Scaling ex- 
ponents are universal: they do not depend on the actual 
value of the forcing and boundary conditions, and they 
are fully characterized by the order of the anisotropy. 
Non-universal effects are felt only in coefficients multi- 
plying the power laws. 

Similar problems, like the existence of scaling laws in 
anisotropic sectors and, if any, the values of the corre- 
sponding scaling exponents are at the forefront of ex- 
perimental, numerical and theoretical research for real 
turbulent flows. Only few indirect experimental investi- 
gations of scaling in different sectors |a, p| and direct de- 
composition in numerical simulations [M-B , 11 1 have been 
attempted up to now. 

The question is still open, evidences of a clear improv- 
ing of scaling laws by isolating the isotropic sector have 
been reported, supporting the idea that the undecom- 
posed correlations are strongly affected by the superpo- 
sition of isotropic and anisotropic fluctuations [^|ll|. On 
a theoretical ground, only recently it has been highlighted 
the potentiality of 5*0(3) decomposition to quantify dif- 



ferent degrees of anisotropics for any correlation function 
Q. On the basis of this analysis, preliminary experimen- 
tal evidences of the existence of a scaling law also in 
sectors with total angular momentum j = 2 have been 
reported |^,^. The value of the exponent for the second 
order correlation function being close to the dimensional 
estimate C^~^(2) = 4/3, |0 (where, from now on, sub- 
script d denotes the dimensional value). 

Tipically, experimental investigations in real turbu- 
lent flows are flawed by the contemporary presence of 
anisotropics and strong non-homogeneities. The mean- 
ing of SO {3) decomposition becomes opaque in presence 
of strong non-homogeneities and also the very existence 
of scaling laws cannot be given for granted [^sj. 
To overcome such difficulties, some of us performed 
(see 0) the numerical investigation of a "Random- 
Kolmogorov Flow" (RKF), a fully periodic Kolmogorov 
flow with random forcing phases, (5-corrclated in time. 

In this Letter, we present a more extended analysis 
of the same data set, but focusing on new evidences 
that anisotropic scaling exponents are indeed universal 
and anomalous, i.e. they do not follow simple dimen- 
sional scaling. In order to do this, we also give a clear 
phenomenological background able to predict the dimen- 
sional scaling in anisotropic sectors. The Letter is orga- 
nized as follows. First, we present a simple dimensional 
argument for all anisotropic sectors of structure functions 
of any order. With respect to this dimensional prediction, 
we show that anisotropic exponents are indeed anoma- 
lous. Moreover, we show that by performing a random re- 
shuffling of the velocity phases (in the stationary ensem- 
ble of anisotropic configurations), the leading anomalous 
scaling is filtered out and the sub-leading dimensional 
prediction is recovered. This is both a test of our dimen- 
sional prediction and a clean indication that forced veloc- 
ity correlations are dominated by inertial terms; in this 
sense one may refer to them as the "equivalent" of the 
zero-modes responsible for anomalous scaling and uni- 
versality in linear hydrodynamical problems ]T3[ |. These 
findings leads to conclude that anisotropic fluctuations 
in turbulence are anomalous and universal. 
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We recall few details on the numerical simulations 0. 
The RKF is fully periodic; the large-scale anisotropic 
random forcing points in one direction, z, has a spatial 
dependency only from the x coordinate and it is different 
from zero at the two wavenumbers: h\ = (1,0, 0),fc2 — 
(2,0,0). Namely, SS^{i,2]) = ^»,3/{i,2} exp (6'{i,2}), 
where /{i,2} are fixed amplitudes and ^'{1,2} are inde- 
pendent random phases, 5-correlated in time. Random 
phases gives an homogeneous statistics, without destroy- 
ing the high anisotropy introduced by the chosen forced 
wavenumber. We simulated the RKF at resolution 256^ 
and collected up to 70 eddy turn over times. 

Anisotropy is studied by means of 5'0(3) decomposi- 
tion of longitudinal structure functions : 



{v(x + R) -v{x))-R 



(1) 



where we have kept only the dependency on R and ne- 
glected the small non-homogeneous fluctuations. We ex- 
pect that the undecomposed structure functions are not 
the real "scaling" bricks of the theory. Theoretical and 
numerical analysis showed [^0,0 that one must first 
decompose the structure functions onto irreducible rep- 
resentations of the rotational group and then study the 
scaling behavior of the projections. In practice, being the 
longitudinal structure functions scalar objects, their de- 
composition reduces to the projections on the spherical 
harmonics: 




10 100 
FIG. 1. Log-log plot of the 4-th order structure function 

projections, 15*4™ (-R)|, in all sectors {j,m) with a high sig- 
nal-to-noise ratio. Notice that the isotropic sector remains 
the most energetic and the clear improving of scaling proper- 
ties if measured on the projected quantities 

dtUi + UkdkUi + UkdkUi + UkdkUi = ~d,p -I- j/Aui. (4) 

The major effect of the large-scale field is given by 
the instantaneous shear Sik = dkUi which acts as an 
anisotropic forcing term on small scales. A simple di- 
mensional reasoning can be done as follows. Let us first 
consider the equation of motion for two point quantities 
{ui{x')ui{x)) in the stationary regime; we may balance 
inertial terms and shear-induced terms as follows: 



Sj,{R)^Yl E S^^- {\R\)Y,^{R) 

j=0 m=~j 



(2) 



As usual, we use indexes (j, m) to label, respectively, the 
total angular momentum and its projection along a ref- 
erence axis, say z. The whole physical information is 
hidden in the functions Sp^{R). In particular, the main 
question we want to address here concerns their scaling 
properties: 



A,^\Rf^Pl 



(3) 



First, we need an estimate for the "dimensional" val- 
ues of the exponents Cd (p) all sectors. Our argument is 
based on the idea that large-scale energy pumping and/or 
boundary conditions are such as to enforce a large-scale, 
anisotropic, driving velocity field U. Dimensional pre- 
dictions for intermediate (small) scales anisotropic fiuc- 
tuations may then be obtained by studying the influence 
of the large-scale U on the inertial range. Let us there- 
fore evaluate the weight of anisotropic contributions as 
it comes out from a balance between inertial advection 
of the small scales and the "shear effect" , induced by the 
instantaneous large scale velocity configuration. Decom- 
posing the velocity field in a small scale component, u, 
and a large scale, strongly anisotropic component, U, one 
finds the following equation for the time evolution of u: 



{ui{x')uk{x)dkUi{x)) ^ {S.ik{x)ui{x')uk{x)) , (5) 

which allows for a dimensional estimate of the anisotropic 
components of the LHS in terms of the RHS shear in- 
tensity and of the {uu) isotropic part. Similarly for 
three point quantities we have (neglecting tensorial de- 
tails): (uuudu) ~ (Suuu) which can be easily general- 
ized to any order velocity correlations. The shear term 
is a large-scale "slow" quantity and therefore, as far as 
scaling properties are concerned, we may safely estimate: 
{Sik{x)ui{x')uk{x)) D^k{ui{x')uk{x)). Here the ma- 
trix Dik is associated to the combined probability to have 
a given shear and a given small scale velocity configura- 
tion. Clearly the Dik tensor brings angular momenta 
only up to j = 2. One may therefore argue, by using 
simple composition of angular momenta, the following 
dimensional matching |L9||: 



SUR) ^ RS 



(6) 



where on top of each term we have written the total an- 
gular momentum of that contribution. 5^(i?) is the pro- 
jection on the j-th sector of the p-th order correlation 
function at scale R (see equation and S is the inten- 
sity of the shear term, Dik, in the j = 2 sector. 
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2. Comparison between the dimensional estimate, 

+ (straight hnes), the measured exponents, 

o) and the exponents, Q.(j)\ obtained after random de- 

(x), for p = 2, 4, 6. Top: sector j = 6, bottom: sector 



For instance, the leading behaviour of the j = 2 
anisotropic sector of the 3-t/i order correlation function 
in the LHS of (^) is given by the coupling between the 
j = 2 components of Z?ife and the j = sector of 
the 2-th order velocity correlation in the RHS of (||): 
SliR) - RS ■ Sl{R) - i?«r'(3). By using the same 
argument and considering that now we know the scal- 
ing of j = and j = 2 sectors of the third order cor- 
relation, we may estimate the scaling exponents of the 
fourth order correlation for j = 2,4. From equation (|^), 
we have the dimensional matching in the j = 2 sector: 
SliR) RS- S^iR) i?'^r'(4) and in the j = 4 sector: 
SiiR) ^ RS ■ SiiR) - R^S^ ■ 55 (i?) - i?«r'(4). The 
procedure is easily extended to higher orders: 



cr 
cr 



ip) 



cr\p- 
cr'ip 
cr'iP- 



l) + l = (p + 2)/3, p>2 
l) + l = (p + 4)/3, p>3 
1) + 1 = (p + 6)/3, p>4 



(7) 
(8) 
(9) 



which can be summarized as 



Cip) = 



(p+j) 



where intermittency effects in the isotropic sector have 
been neglected for simplicity. In this way, giving as in- 
put only the isotropic exponents, Cd~'^ip)' ^'"^ able to 
predict the scaling exponents up to j = 2 for the third 
order structure functions, to j = 4 for the fourth order, 
to j = 6 for the fifth order and so on. We may do a 
little better by giving a prediction also for anisotropic 
fluctuations of second order correlation functions. This 
cannot be simply obtained by using the equations of 
motion, because the first one involving velocity corre- 
lations at different spatial locations, i.e. inertial range 
quantities, is that for dt {ui{x)uj{x')), which fixes a con- 
straint only for the third order correlation function (j^). 



A way out is to ask the second order anisotropic fluc- 
tuations to be analytic in the shear intensity, S, con- 
sistently with what one flnds for higher order structure 
functions by the above dimensional estimate. With this 
assumption, we recover for j = 2 Lumley prediction [ pT[ , 
C^~^(2) — 4/3 by simply writing the first two terms di- 
mensionally consistent with an expansion in the shear 
intensity: {uu) ~ (ei?)^/^ + SR^/^ + ... where the first 
corresponds to the isotropic scaling, while the second cap- 
tures anisotropics up to j = 2 (higher j-sectors could be 
captured by adding other terms in the expansions). By 
using this argument, we may now remove the limit of va- 
lidity of the dimensional prediction, (@-^, and extend it 
to all p values. 

We now come to our numerical results for the S'0(3) de- 
composition of longitudinal structure functions. In Fig- 
ure ^ we present for the A-th order longitudinal structure 
function, an overview of all sectors (j, m) which have a 
signal-to-noise ratio high enough to ensure stable results. 
Sectors with odd js are absent due to the parity symme- 
try of our observable. We measured anisotropic fluctua- 
tions up to J =6. Scaling exponents can be measured in 
almost all sectors except for j = 2 where an annoying os- 
cfllation in the sign of Sl'™'{\R\) prevents us from giving 
a quantitative statement. 

We notice, as it is also summarized in Figure |[ that all 
measured exponents show a clear departure from the di- 
mensional prediction. For example we measure in the j — 
4 sectors the values: (^^(2) = 1.65(5), C'*(4) = 2.20(5), 
C,^{Q) = 2.55(10), and in the j = 6 sector: (^(2) = 3.2(2), 
C^(4) = 3.1(2), C®(6) = 3.3(2). This is a flrst clear sign 
that anisotropic scaling exponents are intermittent. 
The importance of being anomalous does not stand on 
the exact values of the exponents, but on the connection 
between anomalous scaling and universality. Indeed, if 
correlation functions in the inertial range are not given 
by a dimensional matching with the large-scale shear, it 
means that they are fixed only by the inertial part of the 
Navier-Stokes evolution. In other words, they should en- 
joy strong universality properties with respect to changes 
of the large-scale physics, similarly to what happens to 
"zero-modes" responsible for anomalous scaling and uni- 
versality in linear hydrodynamical problems [|l3| . 
Such a statement can even be tested in a different way. 

We have taken the stationary configurations of the 
RKF and randomly re-shuffled all velocity phases : 
Ui{k) Pii{k) ui{k) exp{i 9i{k)), where Pii{k) is the in- 
compressibility projector and 9i{k) = —9i{—k). In this 
way we expect to filter out the dominant intermittent 
fluctuations coming from the inertial evolution, or at 
least those intermittent contributions connected to non- 
trivial phase organization. The rationale of the above 
statement comes from the observation that anomalous 
scaling, in linearly advected hydrodynamical models, is 
connected to the existence of statistically preserved struc- 
tures with highly complex geometrical properties [E0[ . 
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FIG. 3. Comparison of scaling properties before (+) and 
after (x) phases randomization of the 4-f/i order structure 
function for j = 4. Straight lines are the best fit slopes in 
the inertial range. In the inset, the changes for the logarith- 
mic local slopes (same symbols); the horizontal dashed line 
corresponds to the dimensional prediction ^^(4) = 8/3. 

We imagine that once canceled the anomalous fluc- 
tuations, the subdominant fluctuations, due to the di- 
mensional balancing with the forcing-shear terms, should 
show up. Still, it is worth to remark, the statistics of the 
velocity field stays non-gaussian. 

In Figure ^ we show the results for the decomposition 
of A-th order structure functions (after phase random- 
ization) in the j = A anisotropic sector. As it can be 
seen, scaling properties change significantly going from 
the anomalous value (before randomization) to the di- 
mensional predictions (after randomization). 
This happens for all sectors and moments wc have mea- 
sured, as it is summarized in Figure |^, with the notable 
exception of the second order structure function where 
phase randomization has almost no effect. An interesting 
fact which can have two explanations. Phases randomiza- 
tion is not enough to completely filter out intermittency, 
especially for two points quantities which should be less 
sensible to phase correlation. Or, as noticed before, be- 
cause second order correlation function is not constrained 
by any equation of motion, dimensional scaling may never 
exist for it even not as a sub-leading contribution. This 
is an important point which certainly deserves further 
numerical, experimental tests. 

In conclusions we have presented a dimensional ar- 
gument able to predict, by means of a matching be- 
tween inertial correlations and shear-induced inertial 
terms, scaling exponents for all structure functions in 
any anisotropic sector. We have shown by a direct nu- 
merical simulation that anisotropic scaling exponents de- 
viate from the previous dimensional prediction, show- 



ing anomalous values. When performing a random re- 
shuffling of all velocity phases, the dimensional scaling 
comes out as a sub-leading contribution. Everything 
points toward the conclusion that anisotropic fluctua- 
tions are anomalous and universal. 
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